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ABSTRACT 



Context. Strong outbursts in very young and embedded protostars are rare and not yet fully understood. They are believed to originate 
from an increase of the mass accretion rate (M„cc) onto the source. 

Aims. We report the discovery of a strong outburst in a low-mass embedded young stellar object (YSO), namely ^mass-j0542484^ - 
0816347 or [CTF93]216-2, as well as its photometric and spectroscopic follow-up. 

Methods. Using near- to mid-IR photometry and NIR low-resolution spectroscopy, we monitor the outburst, deriving its magnitude, 
duration, as well as the enhanced accretion luminosity and mass accretion rate. 

Results. [CTF93J216-2 increased in brightness by ~4.6, 4.0, 3.8, and I.9mag in the J, H, bands and at 24yum, respectively, 
corresponding to an Ljo; increase of ~20Lq. Its early spectrum, probably taken soon after the outburst, displays a steep almost 
featureless continuum, with strong CO band heads and HtO broad-band absorption features, and Bry line in emission. A later spectrum 
reveals more absorption features, allowing us to estimate reff~3200 K, M,~0.25 Mq, and M„cc~ 1-2x10"'' Mq yr"' . This makes it one 
of the lowest mass YSOs with a strong outburst so far discovered. 
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■ 1. Introduction 

00 Most of the stellar mass in low-mass YSOs is assembled within 
the first 10^ yr of t heir evolution (i. e. class 0, 10** yr, and Class I , 
O ; lO^yr: see e.g., iLada & Wilking [19841 lAndre et alj [19931) . 
p^j • During this stage, the YSO luminosity is thus expected to be 
T— I ' dominated by accretion. However, several st udies, including the 
latest Spitzer Space Telescope surveys (e. g..l Enoch et al.|[2009t 
[Evans et al. 2009), have found that more than 50% of the embed- 
^ ^ ded YSOs have Lhoi and Mace values considerably lower than 
• I— I , those theore tically predicted (i.e. ~2xlO~ ^MHyr~' for solar- 
X ■ mass YSOs: [Shul[l977t iTerebev et alJll984 and roughly of the 
JH ' same order of magnitude as the Classical T Tauri s tars (CTTs; 
' i.e. 1 0'*^ < Mace < 10 ^' vr ' : e. s.. IWhiteXH illenbrand 



l2004l) . Among several hypotheses, a likely explanation is that 
the mass accretion is episodic, and the protostars with the 
lowest luminosities are t h ose observed in quiescent accretio n 
states (lEnoch et all 120091: lEvans et al.l 120091: IVorobvovl l2009h . 
Non-steady mass accretion is often observed in CTTs, such as 
EXors and FUors (lasting from a few months to several decades), 
in which Mace incre ases by several orders of m agnitude up to 
M„cf~10 "^Moyr"^ ( iHartmann & Kenyonlll996l) . It is thus rea- 



* Based on observations collected at the ESO/NTT (082.C-0264), at 
the REM telescope La Silla, Chile, and at the the Italian Telescopio 
Nazionale Galileo (TNG), operated on the island of La Palma by 
the Fundacion Galileo Galilei of the INAF (Istituto Nazionale di 
Astrofisica) . 



sonable to believe that a similar mechanism also exists in earlier 
and more embedded YSOs. Unfortunately, there is little direct 
observational evidence of outbursts in Class I YSOs, and so far, 
only a few clear cases have been detected (e. g. V 1 647 0ri out- 
bursts in 2003 and 2008 o r OO Ser in 1995: see e. g. lFedele et all 
l2007HK6sDar et al.l[2007l) . To reconcile theory with observations 
and improve the quality of the statistics, it is mandatory to detect 
and study these rare events. 

With this aim, we started a long-term project to monitor 
the NIR flux and spectroscopic variability of embedded YSOs 
(mostly Class I and Flat sources) in nearby, young, and ac- 
tive star-forming regions (namely L 1641, CrA, and the Serpens 
Molecular Cloud). This letter reports the outburst of an embed- 
ded YSO in L 1641, namely 2mass-j0542484B-0816347, here- 
after [CTF93]216-2 (q'2ooo=05"42"'48.48', 5200()=-08°16'34'.'7). 
This object was identified by our group as a low-mass embed- 
ded YSO (spectral type M, circumstellar Av~18mag) with a 
flat spectral index (a'=0.25, derived by fitting all the photomet- 
ric data points from 2.2 to 24 //m) and a bolometric luminos- 
ity of ~1.9Lq (Caratti o Garatti et al. in prep., hereafter CoG). 
It has been nam ed [CTF931216-2 , because it is relatively c lose 
to [CTF93]216 dChen et alJll99l IChen & Tokunagalll994 . lo- 
cated about 38'.'5 SW. Our Spitzer/IRAC images indicate that 
both YSOs have precessing jets, thus they could be part of a wide 
binary system (~17 300AU, assuming a distance (i=450pc). 
During our recent survey in October 2010 with the robotic tele- 
scope REM (see Section|2]i, we detected for [CTF93]216-2 a 
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brightness increase of several magnitudes with respect to the 
2MASS images. We then compared our new images with the 
acquisition image in the band and the NIR spectrum of this 
source acquired in February 2009, discovering that the outburst 
was already in progress. 



2. Observations 

NIR spectra were obtained at the ESO-NTT with 
SofI dMoorwood et al.1 IT998I) (on the 13 Feb. 2009) and 
at the 3.5-m Italian te lescope TNG (13 Oct. 2010) with 
NICS jBaffa et al.l 1200 Ih . adopting the usual ABBA config- 
uration. The SofI spectrum was taken with the red grism 
(1.51-2.5yum), a a.'6 sUt H -1000), and a total integration 
time of 1800 s. The full width half maximum {FWHM) in the 
dispersion direction, measured from Gaussian fits to the OH 
sky lines, was ~19A (-260 kms ') in the K band. The NICS 
low -resolution spectrum was acquired using a slit width of 1" 
(•R -600) for the JH and HK grisms (1.15-1.75jum and 1.4- 
2.5 //m, respectively) with a total integration time {Int) of 480 
and 160 s for the JH and HK grisms, respectively. The measured 
FWHM in the dispersion direction was -33 A (-600kms"') 
in the K band. Telluric and spectrophotometric standards 
were observed to correct for the atmospheric transmission and 
flux-calibrate the spectra. The wavelength calibrations were 
performed using a xenon lamp in the infrared. 

Most of our J, H, and Ks images were collected between 
the 10 October a nd the 6 Novemb er 2010 at the 60-cm robotic 
telescope REM dZerbi et alJl200lL ESO La Silla observatory) 
with the NIR camera REMIR jConconi et al .1 12004') and a 150 s 
total integration time per filter. An additional H band image 
was obtained with NICS (13 Oct. 2010) with a 18 s total inte- 
gration time. Additional K s photometry (5 Nov. 2003) with the 
near-IR instrument UIST dRamsav Howat et al.l l2004l) was re- 
trieved from the UKIRT data archive (Int-15 s). Moreover, an- 
other photometric data point was derived from the K, band of 
the SofI Kg acquisition image {Int-6 s). Early epoch photometry 
was retrieved from th e Two-Micron All Sky Survey (2MASS; 
ISkrutskie et al] l2006l J, H, K,, band; November 1998) and 
from the DEep Nea r-Infrared southern sky Survey (DENIS; 
lEpchtein et alii 19971 J and Ks bands; January 1997). All the 
raw data were reduced using IRAF packages, applying standard 
procedures for sky subtraction, dome flat-fielding, and bad pixel 
and cosmic ray removal. Photometric calibration was obtained 
by means of photometric standard stars, except for the SofI Ks 
acquisition image and the UIST image, where 2MASS photom- 
etry of field stars was used. 

Finally, we also used additional Spitzer/IRAC photometric 
data (3.6, 4.5, 5.8, and 8 yum, obtained on the 8 Oct. 2004), MIPS- 
24jum (taken on the 2 Apr. 2005 and 27 Nov. 2008), MIPS 70 
and I60//m (2 Apr 2005), as well as Spitzer/IRS low -resolution 
spectroscopy (5.2-39yum, obtained on the 12 Mar. 2007), pre- 
sented in CoG. Additional SCUBA/JCMT pho t ometry (at 450 
and 850 jum) was taken from lDi Francesco et al. I (l2008l) . 



3. Results 

3.1. Photometry 

In Figure [1] pre- and outburst images of [CTF93]216-2 in the 
Ks band (top panels) and the Spitzer/MIPS 24 //m band (bot- 
tom panels) are shown, clearly indicating an increase in the 
object brightness. Lightcurves in the J, H, Ks bands, and at 
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Fig. 1. Ks (top panels) and MIPS -24 /^m (bottom panels) im- 
ages of [CTF93]216-2 and its surroundings before (left) and af- 
ter (right) the outburst. The position of the YSO [CTF93]216 is 
also indicated. 

Table 1. Photometry of [CTF93]216-2. 



Date 


JD 


J 


H 


Ks 


24 ^im 


(d.m.y) 


(2450000+) 


(mag) 


(mag) 


(mag) 


(mag) 


10.1.1997 


0458.6 


16.1±0.2 




I1.55±0.I 




18.11.1998 


1135.7 


16.4±0.1 


13.53±0.02 


11.76±0.03 




5.11.2003 


2948.8 






10.9±0.2 




2.4.2005 


3462.8 








2.48±0.1 


12.3.2007 


4172.0 








2.49+0.06° 


27.11.2008 


4798.2 








0.58+0.04 


13.2.2009 


4875.5 




9.55±0.06'' 


7.98±0.05 




10.10.2010 


5479.9 


11.77±0.08 


9.66±0.05 


8.13±0.07 




12.10.2010 


5481.9 


11.83±0.09 


9.64±0.04 


8.12±0.05 




13.10.2010 


5482.7 




9.70±0.05 






19.10.2010 


5488.9 


11.82±0.07 


9.63±0.06 


8.15±0.07 




28.10.2010 


5497.8 


12.01 ±0.1 


9.76±0.09 


8.25±0.1 




06.11.2010 


5506.6 


12.00±0.07 


9.70±0.06 


8.24±0.08 





Notes. Photometry derived from the Spitzer-IRS flux-calibrated 
spectrum. Photometry derived from the SofI flux-calibrated spec- 
trum. 



24 yum are shown in Figure |2] and the corresponding photo- 
metric data points along with their uncertainties are reported 
in Table [T] We stress that the photometric data points at 24 jim 
(on the 12 Mar. 2007) and in the H band (on the 13 Feb. 2009) 
were derived from the Spitzer/IRS and from the NTT/SofI flux- 
calibrated spectra, respectively. The 24jum data points (Fig. |2] 
upper panel) give us some constraints on the date of the outburst, 
i. e. after March 2007 and before November 2008. Pre-outburst 
photometry from DENIS, 2MASS, and the UIST Ks images in- 
dicates that the object is variable (in both J and Ks bands, we 
have just one 2MA SS data point for the H band), as is often 
the case for YSOs dCarpenter et al.ll200lh . Comparing 2MASS 
(November 1998) with our earliest outburst SofI photometry 
(February 2009), we measure the amplitudes A/rj-3.8 mag, and 
A//-4mag, corresponding to a flux increase AF by a fac- 
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Fig. 2. [CTF93]216-2 lightcurves. J, H, K^, and 24 /jm data are 
represented as black squares, red dots, blue triangles, and green 
diamonds, respectively. 
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tor of -33 and 40, respectively. Considering our earliest post- 
burst photometric point in the J band (October 2010), we de- 
rive A7~4.6mag (or AF~69), whereas the Amag at 24 //m is 
~ 1.9 mag (AF~6). Our data do not allow us to determine the 
exact date of the peak, but it is clear that the NIR lightcurves are 
still close to the maximum, and that this plateau phase has lasted 
at least 2 years, with a small decrement in the H and K, bands 
of about 0.1 mag between February 2009 and October 2010. 
Moreover, our latest photometry (November 2010) seems to in- 
dicate that the rate of decrease, between October and November 
2010, has increased to ~0. 2 mag/month in the J band and 
~0.1 mag/month in the H and K, bands. Pre- and outburst spec- 
tral energy distributions (SEDs) are shown in Figure [3] where 
all available photometric and spectroscopic data (from this work 
and from CoG) are reported. We note that both the spectral in- 
dex and bolometric luminosity have change d, as already reporte d 
in other YSO outbursts (see e.g. OO Ser. iKospal et al.ll2007l) . 
During the outburst, the YSO SED has become bluer and flat- 
ter, and the a value, computed between 2.2 and 24 /im, var- 
ied from 0.69 to -0.04. To compute Lboi, and estimate the out- 
burst Amag in the Spitzer/IRAC photometry, we fit the dered- 
dened Amag amplitude from 1.25 to 24 /^m by means of a power 
law, obtaining Am(3.6/im)=2.9mag, Am(4.5jum)=2.6mag, 
A;7i(5.8/im)=2.3 mag, and AOT(8yum)=2.1 mag. We also assume 
that the SED flux experienced no significant increase for 
A >24yum. As a result, we estimate an Li,oi value of ~22Lq 
during the outburst, i. e., ALboi ~20Lo with respect to the pre- 
outburst state. 

3.2. Spectroscopy 

Our Sofl and NICS low resolution spectra (taken in February 
2009 and October 2010, respectively) are shown in Figure |4] (in 
black and red). 

The Sofl spectrum shows an almost featureless continuum, 
with strong veiling and IR excess. The absorption water band 
between 1.7 and 2.1 fim, typical of M spectral types with a 
Teff <3800K, is clearly visible. In addition, deep CO band- 
head absorption lines are detected at 2.29/im (2-0), 2.32//m (3- 
1), and 2.35 yKm (4-2). Finally, Bry emission, usually associated 
with accretion, is observed in our Sofl spectrum. The measured 



Fig. 3. Spectral energy distribution of [CTF93]216-2. Pre- 
outburst data: black triangles (photometry), and black line 
(Spitzer-IRS mid-lR spectrum). Outburst data: red squares (pho- 
tometry), and red line (Sofl NIR spectrum). 

Bry flux is 9.6(±1.8)xl0-'4erg s"' cm^^ (EW~-3.6A), peaking 
at 2.166yum. 

The flux of the NICS spectrum, taken 607 days later, shows 
a slight decrease. The continuum still has a similar shape, but 
the absorption H2O band depth is more prominent. The EW of 
the CO fines slightly increases by ~1 A (9.8 to 10.6 A and 7.2 
to 8.2 A, for the 2-0 and 3-1 band heads, respectively). More 
CO lines at 2.38 /urn (5-3), and 2.41 /zm (6-4) are detected. In 
addition, a few more narrow-band absorption features become 
visible on the continuum, i.e. Cai (2.26 /^m), the Nai doublet 
(2.20-2.21jum), Mg i (1.59- 1.71 //m), and the Al i doublet (1.67 - 
1.68 fim). No relevant features are detected in the J band, except 
for the absorption H2O band between 1.3 and 1.5jum. The Bry 
line is not resolved, peaking at 2.166yum. The measured flux is 
9(+2)xl0"''*ergs"' cm"- (EW~-3.7A). We note that no signa- 
ture of any emission line (such as H2 or [Fe 11]) from the jet has 
been detected on-source in both spectra. 

To more tightly constrain the YSO spectral type, we sepa- 
rately compared limited regions of our spectra (in the H and K 
bands) to model spectra. Veiling and scattering (not taken into 
account by our fit) are expected to smoothly and slowly vary 
as a function of wavelength and thus, they shou ld not affect 
the sp ectrum in individual bands significantly (e. g. JScholz et al.l 
l2010h . For our comparison, we us ed a series of AMES-DUSTY 
model spectra dAllard et al.l 1200 ih with Teg ranging from 2500 
to 3900 K and logg=4.0, as expected for young stars. We var- 
ied Teff to obtain a consistent solution for the two bands, using 
the Ai/=18 mag derived by CoG. We find a reasonable match be- 
tween observed and model spectra for reff=3200±200K. As an 
illustration, the model for 3200 K is shown in Fig.|4](blue spec- 
trum), shifted by -0.3x10"'" ergs"' cm"-/im"' for clarity. 

4. Discussion 

Our discovery of the [CTF93]216-2 outburst gives us a rare op- 
portunity to study boosted accretion in young embedded proto- 
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1.2 1.4 1.6 1.8 2.0 2.2 2.4 



Wavelength (/itm) 

Fig. 4. Flux calibrated MR spectra of [CTF93]216-2: SofI spectrum (black), and NICS spectrum (red). F or comparison, we show 
an AMES-DUSTY model spectrum for effective temperature of 3200 K and log g=4.0 dAllard et alJl200lh . The DUSTY spectrum 
has been shifted by -0.3x 10 '° erg s ' cm"^ yum ' for clarity. The detected features are also labelled. 



stars, probing whether episodic mass accretion can reconcile the 
low accretion rates observed in young embedded protostars with 
the theoretical predictions. 

The Bry luminosity, corrected for the circumstellar ex- 
tinction, is often used to derive an estimate o f the accre- 
tion luminosity (see e. g. iMuzerolle et al.l Il998t iNatta et aP 
l2006h . Assuming that Av=18+3mag a nd J=450pc, we derive 
L(Br7)~3.5(0.9)xlO"^^Lo and, from the IMuzerolle et al.l (Il998h 
relationship, we obtain Lflce~22+8 Lq, which is close to the de- 
rived ALfe„/~20LQ and indicates that L^cc ~ t^l^boi- 

On the ot her hand, M^cc can be derived from Mace — Lace x 
l.25R,/GM, ( iGullbring et al.l ll 9981) . where M» and /?» are the 
stellar mass and radius, and the number 1 .25 is derived by as- 
suming a value of 5R^ for the inner radius of the accretion disk. 

Assuming tha t Teif =3200+200 K and L.~1.9+0.1 L©, from 
ISiess et al.1 (|2000|) models, we obtain M. =0.24+0.04 M© and 
Rt=4A±QARQ, which gives a mass accretion rate of 
1.3xlO"'^M0yr"'. While M, is well constrained by Tes, the ra- 
dius depends on L«. Thus this last should be considered as an 
upper limit, since pre-outburst Lace is unknown. However, even 
a 50% decrease in L, would affect R, and thus Mace of just 20%. 

Compared to typical accretion rates of Class I YSOs with 
similar masses (i.e. 10~^M(Tiy r~'; e.g., CoG; IScholz et al.l 
l20Tol:IWhite & Hill enbranJ2004'), the derived Mace is about two 
orders of magnitude higher This value is about an order of mag- 
nitude lower than what some episodic-accretion mod els predict 
for th ese early YSO bursts (i.e lO^^MQyr"'; e. s. IVorobvovl 
I2009h . which should resemble the more evolved FUor counter- 
parts. On the other hand, we note that the magnitude of lumi- 
nosity change during this outburst (~10) is similar to those of 
V 1647 Ori and OO Ser (~8). Our object would probably differ 
from the others in its pre-outburst Lhoi, which is about two times 
lower, because of its l ower mass. Moreover, as already noted 
in previous cases (e. g. iHodapp et al.l 1 19961 iFedele et"aDl2007t 
iKospal et al.1l2007l) . the [CTF93]216-2 outburst shows similari- 
ties with both FUor and EXor events, i. e. a NIR featureless FUor 
spectrum with strong absorption CO band heads, but with Bry 
line in emission as in EXors. The event duration (if confirmed 
by additional observations) and amplitude are in-bet ween those 
of EXors and FUors, thus, as noted by other authors dGibb et al.l 



l2006HFedele et al.ll200'n) EXors and FUors might not be distinct 
categories of eruptive events, but instead part of a continuum of 
outburst events. 
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